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Localized heating of a gas by intense laser pulses leads to interesting acoustic, hydrodynamic 
and optical effects with numerous applications in science and technology, including controlled wave 
guiding and remote atmosphere sensing. Rotational excitation of molecules can serve as the energy 
source for raising the gas temperature. Here, we study the dynamics of energy transfer from the 
molecular rotation to heat. By optically imaging a cloud of molecular superrotors, created with 
an optical centrifuge, we experimentally identify two separate and qualitatively different stages of 
its evolution. The first non-equilibrium “gyroscopic” stage is characterized by the modified optical 
properties of the centrifuged gas - its refractive index and optical birefringence, owing to the ultrafast 
directional molecular rotation, which survives tens of collisions. The loss of rotational directionality 
is found to overlap with the release of rotational energy to heat, which triggers the second stage of 
thermal expansion. The crossover between anisotropic rotational and isotropic thermal regimes is 
in agreement with recent theoretical predictions and our hydrodynamic calculations. 
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Intense ultrashort laser pulses can deposit energy in a 
transparent gas medium either by ionizing [TJ [2] or spin¬ 
ning up[3l [4] its molecular constituents. Heating of the 
gas is followed by the formation of an acoustic wave [5]- 
[8] which leaves behind a long-lived low-density depres¬ 
sion channel [9]. Laser-induced sound emission is key to 
photoacoustic spectroscopy [TOHE], whereas laser control 
of the gas density proved valuable for generating [131 [14] 
and even controlling [4] wave guiding channels in ambient 
air. Because ionization-free delivery of energy by means 
of rotational excitation does not inhibit coherent light 
propagation through the heated gas, it has the potential 
of producing longer wave guides, e.g. if executed inside 
ionization-free filaments [15]. 

Together with the importance of rotational heating, re¬ 
cent developments in controlling molecular rotation with 
nonresonant laser pulses [T6l [TT] have stimulated active 
research on the exchange of energy between a rotating 
molecule and its environment. Transient molecular align¬ 
ment has been first proposed [18] and later implemented 
as a powerful probe of collisional relaxation in a series 
of pioneering experiments [TM2T] . The three relaxation 
steps, associated with (i) rotational decoherence, (ii) ro¬ 
tational reorientation and (iii) rotation-translation (RT) 
thermalization, have been identified in the theoretical 
model [22], yet found to overlap in time too closely to 
enable an individual experimental study of each process 
separately [21]. 

Extreme rotational excitation of molecular 
superrotors [23] with an optical centrifuge [24] [25] 
allowed us to separate the time scales of the relaxation 
mechanisms and investigate them in great detail. In our 
earlier work[26l [27] . we have focused on the rotational 
decoherence of superrotors which takes a few hundred 
picoseconds under ambient conditions and ends well 


before the onset of the rotational reorientation. Ul- 
trahigh values of the molecular angular momentum J, 
provided by the centrifuge, together with the propensity 
of collisions to conserve its orientation [22] results in 
the “gyroscopic stage” [28] , which outlives the rotational 
coherence by a few nanoseconds. Here, we study the 
dynamics of this relaxation step by tracking an optical 
birefringence of the centrifuged gas in the direction 
perpendicular to the direction of the centrifuge. This 
transverse birefringence stems from a strong permanent 
confinement of molecular superrotors in the plane of 
their rotation. In the longitudinal direction along the 
centrifuge, the change in refractive index owing to the 
fast unidirectional molecular rotation creates a refractive 
“gyroscopic channel”, also detected in this work by 
means of the phase contrast imaging. 

As the superrotors undergo collisional relaxation, their 
rotational energy is released, driving the temperature of 
the gas up. The rotation-translation energy transfer is 
accompanied by a crossover from the non-equilibrium 
gyroscopic to the thermal phase of molecular dynamics, 
which to the best of our knowledge is observed here for 
the first time. Hydrodynamics of a locally heated gas re¬ 
sult in the corresponding change in its pressure, followed 
by the emission of an acoustic wave and the formation 
of a low-density thermal channel. By directly measur¬ 
ing the evolution of the gas density on the time scale 
of tens and hundreds of nanoseconds, following its ther¬ 
mal expansion for up to 3 microseconds, we are able to 
project it back to the origin of the thermal phase and 
analyze the heating rate at the RT crossover. The anal¬ 
ysis, based on our full hydrodynamic calculations, points 
at the close temporal overlap between the rotational re¬ 
orientation and rotational cooling, in accord with recent 
theoretical predictions [28]. 
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FIG. 1: Schematic diagram of the experimental setup. An ensemble of molecular superrotors is created in a gas cell with an 
optical centrifuge (red beam). Nanosecond or picosecond probe pulses (green and blue beams, respectively) are delayed with 
respect to the centrifuge pulses and propagate either collinear with, or perpendicular to the centrifuge, creating an image of 
the rotationally excited volume of gas on a CCD camera either in the longitudinal or transverse direction, respectively. DM: 
dichroic mirror, FL: focusing lens, IL: imaging lens, KE: knife edge (shown along y instead of z axis for clarity), LP: linear 
polarizers at ±45° to y axis, F: frequency filter. An example of the longitudinal image with a circularly expanding sound wave 
is shown in the upper right corner. Images in the lower row were taken in the transverse geometry at early (left), intermediate 
(middle) and late (right) time moments. The leftmost image was recorded with the two linear crossed polarizers in place. The 
rightmost schlieren image was recorded with the knife edge in place. 


Our experimental configuration is depicted in Fig|^ 
The details of our optical centrifuge have been discussed 
in an earlier publication [27]. Briefly, the centrifuge 
shaper, built according to the original recipe of Karcz- 
marek et al [24], is followed by a home built TiiSapphire 
multi-pass amplifier boosting the pulse energy to 30 mJ. 
Centrifuge pulses are about 100 ps long and their linear 
polarization undergoes an accelerated rotation, reaching 
the angular frequency of 10 THz by the end of the pulse. 
The centrifuge beam is focused by a weak 1 m focal length 
lens down to a spot size of 90 /im diameter (full width 
at half maximum) inside a gas cell filled with 0.9 bar of 
oxygen. This focusing results in the peak intensity of 
2.2 X 10^^ W/cm^, significantly lower than the ionization 
threshold of O 2 . We monitor the performance of the cen¬ 
trifuge by means of coherent Raman spectroscopy of the 
created molecular superrotors. The Raman setup is de¬ 
scribed elsewhere [23] and not shown in Fig0 for clarity. 

We image the volume of the centrifuged gas both on 
short and long time scales, r < 23 ns and r < 2 /is, re¬ 
spectively. For a short-time scan, picosecond pulses are 
extracted from the same TiiSapphire ultrafast system as 
the centrifuge, and delayed with femtosecond precision 
by means of computer controlled motorized translation 


stages. Long-time scans are executed with nanosecond 
pulses from a separate YAG laser, while their delay is 
controlled electronically. In both cases, single-shot im¬ 
ages are recorded with a CCD camera and typically av¬ 
eraged over 400 laser pulses. 

Two imaging geometries are implemented: longitudinal 
and transverse, with probe pulses propagating collinear 
with, or at 90 degrees to the direction of the centrifuge. A 
sample longitudinal image, taken 2 jas after the centrifuge 
pulse, is shown in the upper right corner of Figj^ The lo¬ 
cation of the central bright spot, surrounded by multiple 
circular interference fringes, corresponds to the position 
of the centrifuge beam. Because of the centimeter-long 
confocal parameter, the fringe pattern stems from the 
diffraction of probe pulses inside the long density depres¬ 
sion channel and does not reflect the true distribution of 
light intensity (e.g. wave guiding) or refractive index [29]. 
Yet owing to the linearity of the weak probe propagation 
and its extended length, the fringe contrast serves as a 
sensitive indicator of the centrifuge-induced changes in 
the refractive index of the gas. An. 

We analyze the time dependence of An by recording 
a series of longitudinal images while scanning the delay 
time between the centrifuge and probe pulses. An angle- 
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FIG. 2: Radial cross-section of images recorded along the cen¬ 
trifuge beam (as in Fig[^b)) as a function of the delay time 
between the centrifuge and probe pulses. Long-time scan in 
panel (a), carried out with ns probe pulses, extends to 3 /iS 
and shows an acoustic wave traveling with the speed of sound 
away from the central core (see Movie 1 in Supplementary 
Material). In panel (b), we zoom in on the first 23 nanosec¬ 
onds (by switching to ps probe pulses) and find a crossover 
between the rotation-induced and thermal channels around 7 
ns, when the image contrast vanishes almost entirely. Sample 
images of both refractive channels are shown in the square 
insets. The crossover is quantified by calculating the change 
in the image contrast and plotted in panel (c). The red line 
on top of the experimental data points is shown to guide the 
eye. 

averaged radial cross-section of the recorded interference 
pattern is calculated and plotted as a function of the 
delay in Figj^a). The plot shows a permanent central 
core and an outgoing sound wave, similar to the recently 
observed dynamics of plasma filaments pMl [29]. Unlike 
the case of plasma, however, zooming in on the first few 
nanoseconds of the filament formation reveals unique de¬ 
tails of its early history. The instantaneous creation of a 
refractive channel by the centrifuge, as seen in panel (b) 
of Figj^ reflects its non-thermal origin. Indeed, we de¬ 
tect the signal already in the first picoseconds after the 
arrival of a centrifuge pulse - too short a time scale to 


allow collisional thermalization. Instead, the collisions 
seem to suppress, rather than enhance An, which disap¬ 
pears almost completely at around 7 ns. From that mo¬ 
ment on, the channel grows back and eventually emits 
the sound wave seen in Figj^a), in accord with the ex¬ 
pected thermal dynamics. A clear crossover between the 
two regions is illustrated in panel (c) of FigJ^ where we 
plot the contrast of the recorded longitudinal images as a 
function of time. To identify the nature of the refractive 
channels before and after the crossover, we switch to the 
transverse imaging geometry. 

Three typical transverse images, recorded on three 
qualitatively different time scales, are presented at the 
bottom of Fig{^ The schlieren technique [30] is employed 
by inserting a knife edge in the focal plane of the imag¬ 
ing lens, to determine the true distribution of the re¬ 
fractive index (more exactly, its d/dy derivative) at long 
delay times. The side view of a narrow channel is free 
of interference fringes, seen in the longitudinal geometry, 
which greatly simplifies the interpretation of the picture. 
Schlieren image (e) on the right shows two high-density 
waves captured 400 /im above and below the low-density 
depression channel in the center. The snapshot (d) in the 
middle was taken close to the very origin of the thermal 
channel around r ~ 20 ns, before the beginning of the hy¬ 
drodynamic expansion of the gas. At even shorter delay 
times of order of, and below, 5 nanoseconds, our phase- 
contrast imaging approaches its signal-to-noise limit. At 
these early moments, we re-gain sensitivity by taking the 
images between two crossed linear polarizers, probing the 
centrifuge-induced birefringence of the gas (leftmost pic¬ 
ture (c) at the bottom of FigQ. 

The difference between the two processes affecting op¬ 
tical properties of the centrifuged gas before and after 
the observed crossover becomes apparent when we com¬ 
pare the time dependence of the detected optical birefrin¬ 
gence with that of the phase contrast. The two observ¬ 
ables, extracted from the images taken with and with¬ 
out the crossed polarizers, respectively, are plotted in 
Figj^ The birefringence is expressed in arbitrary units 
with 1 being equivalent to the polarization rotation an¬ 
gle of about 1 mrad, while the phase contrast is quanti¬ 
fied as the image intensity contrast. From these results, 
the disappearing channel to the left of the crossover in 
Figj^b) can be clearly correlated with the decaying bire¬ 
fringence signal, and therefore associated with an optical 
anisotropy induced by the directional molecular rotation. 
As the molecular superrotors lose their rotational energy 
to heat, this gyroscopic channel dies off while thermal 
effects take over, causing an isotropic change in the gas 
density reflected in the phase contrast images at r > 7 
ns. To make sure that superrotors are responsible for 
the observed optical changes, we repeat the experiments 
with the centrifuge pulses of the same peak intensity but 
lower terminal rotation frequency. Two points inside the 
dashed rectangle in Fig|^ upper square and lower dia- 
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FIG. 3: Birefringence signal (red triangles) and phase contrast 
(blue circles) of the centrifuged gas, retrieved from the trans¬ 
verse images taken with and without the crossed polarizers 
(left and right insets, respectively) and plotted as a function 
of the delay time between the centrifuge and probe pulses. 
Optical birefringence is induced by the centrifuge instanta¬ 
neously and falls off exponentially with a time constant of 3.4 
ns (solid black curve, fit at r > 100 ps, i.e. after the end 
of the centrifuge pulse). A crossover to the growing phase 
contrast is observed around 7 ns, similarly to the crossover 
in FiglD Two data points in the lower right corner (within 
the dashed rectangle) indicate the drop in the phase contrast 
for the slower rotating centrifuge. The dashed black line is a 
fit to the hydrodynamic numerical calculations of An(T) with 
an RT exchange rate of 3.4 ns, whereas the dashed magenta 
and green lines represent the best fits with that rate constant 
being five times shorter and longer, respectively. The hori¬ 
zontal dash-dotted line indicates the noise floor for the phase 
contrast measurement. 


mond, correspond to the centrifuge producing only slow 
rotors and no rotors at all, respectively. 

To verify the mechanism behind the rotation- 
translation crossover, we numerically simulate the hydro¬ 
dynamics of the ideal gas exposed to a known heat source 
(see Methods). The hydrodynamics equations, corre¬ 
sponding to the conservation of mass, momentum and 
energy in the compressible gas flow, are solved for the gas 
density p, velocity v and pressure P, assuming cylindrical 
symmetry and equilibrium initial and boundary condi¬ 
tions. Rotation-translation energy exchange is modeled 
by adding an external pressure source, ^ {dP/dt)-^rj. = 
1^0 the corresponding differential equation 
for P. As we argue in Methods, an exponential increase 
of the gas temperature (and hence, its pressure in the first 
few nanoseconds) with the same time constant = 3.4 
ns as the decay of the birefringence signal, stems from 
the proportionality of both the birefringence signal and 
the rotational energy to In the expression above. 
To (Pq) is the ambient temperature (pressure) and AT is 
the temperature increase at the end of the thermalization 
process. To estimate AT, we note that the peak intensity 


of our centrifuge pulses is sufficient to spin adiabatically 
only 2% of O 2 molecules occupying the lowest rotational 
state at room temperature. After these molecules are 
excited to J ~ 91 by the centrifuge (as determined by 
means of Raman spectroscopy[23]), each oxygen super¬ 
rotor carries 1.5 eV of rotational energy. Redistributing 
this energy among the whole molecular ensemble and all 
5 degrees of freedom results in AT = 140 K. 

The main result of our calculations is shown in 
Figg a), where the change in the gas density, Ap(p,r), 
is plotted as a function of both the distance from the 
centrifuge and the time since its arrival. One can see 
the density depression channel forming in the center in 
the first 200 ns and an acoustic density wave spreading 
radially away from it. Given the linear dependence of 
the refractive index of a gas on its density, we compare 
the experimentally observed phase contrast with the nu¬ 
merically calculated \Ap{y = 0,r)| in Fig(blue circles 
and dashed black line, respectively). The only fitting pa¬ 
rameter is the proportionality coefficient between the gas 
density and the image contrast. An excellent agreement 
between the experimental and numerical line shapes sup¬ 
ports our interpretation of the gyroscopic channel and its 
role as the energy source for the thermal gas expansion. 
If the RT exchange rate is assumed to be very different 
from the birefringence decay rate r^, the numerical re¬ 
sults do not fit the experimental data as well (dashed 
magenta and green lines in FigJ^. 

Further comparison between our experimental findings 
and numerical model is presented in Figj^ In panel 
(b), the schlieren signal at 1.2 ps is calculated by tak¬ 
ing the derivative of Ap{y, r = 1.2 ps) with respect to 
and plotted together with the y cross-section of the ex¬ 
perimental schlieren image taken at that delay time (cf. 
FigQe)). The dependence of schlieren signals on both 
y and r is shown in FigJ^c) and (d) for the experimen¬ 
tal and numerical data, respectively. The experimental 
scan, carried out with ns probe pulses, extends to 2 ps 
and shows an acoustic wave traveling with the speed of 
sound away from the central core (see Movie 2 in Supple¬ 
mentary Material). The equivalent numerical scan dis¬ 
plays very similar behavior, supporting the validity of the 
model. 

To summarize, we report on the experimental obser¬ 
vation of a refractive “gyroscopic channel”, created by 
an ultrahigh rotational excitation of molecules with an 
optical centrifuge. Owing to the non-thermal nature of 
the channel, it affects the optical properties of a gas, 
such as its index of refraction and birefringence, on ul¬ 
trashort (picosecond) time scale which may prove useful 
for ultrafast optical switching. The gyroscopic channel 
is found to serve as a precursor to the formation of a 
conventional thermal channel (typically originating from 
the laser-induced plasma at much higher laser intensi¬ 
ties), and a supplier of energy to the latter. A crossover 
between the rotation relaxation and thermal expansion. 
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FIG. 4: Comparison of the experimental results with the numerical calculations of gas hydrodynamics. The calculated change 
in the gas density, Ap, is shown in panel (a) as a function of time and distance. Experimentally determined parameters of 
the observed gyroscopic channel are used to simulate the heat source, which initiates the dynamics at point (y = 0,t = 0). 
In panel (b), the derivative of the calculated density prohle d/dy[p(y)] at 1.2 ps is compared with the y cross-section of the 
schlieren image recorded 1.2 ps after the centrifuge (dashed red and solid blue curves, respectively). The dependence of the 
measured and calculated schlieren signals on both space and time is shown in panels (c) and (d), respectively, with the two 
lines indicating the two cross-sections displayed in plot (b). 

driven by the rotation-translation energy exchange, is 
discovered experimentally and reproduced in the numer¬ 
ical gas hydrodynamics calculations. 


Providing a few electron-Volts of rotational energy 
per molecular superrotor, comparable to the energy of 
free electrons in plasma filaments, an optical centrifuge 
is shown to be a non-ionizing alternative for initiating 
and studying thermal processes in ambient gases, e.g. 
photo-acoustics [To] or wave guiding [131 (HI? especially 
valuable because of its capability to precisely control the 
amount of energy transferred to the gas [23] . Translated 
to the control of the optical properties of gases, demon¬ 
strated here, this may be of use in applications ranging 
from ultrafast polarization switching [3T1 - I33] to remote at¬ 
mosphere sensing [34] . Ultrafast rotation of centrifuged 
molecules results in the increasing adiabaticity of molec¬ 
ular collisions and correspondingly slower release of ro¬ 
tational energy. A series of fascinating phenomena, such 
as an explosive thermalization and anisotropic diffusion, 
have recently been predicted for a gas with high concen¬ 
tration of molecular superrotors [28]. This work is a step 
towards exploring this new area of gas dynamics. 


METHODS 

In this work, we numerically solve the following set of 
hydrodynamic equations [35]: 

d/dt{p) = -V[pv], 
d/dt{v) = -V[pv + P], 
d/dt{£) = -V[fv + Pv + q], 

where p, p = pv and £ = \pv‘^ + 5 are the mass, mo¬ 
mentum and energy of the unit gas volume; v its velocity 
and 5 = pCvT its internal energy, with being the heat 
capacity at constant volume and T the gas temperature. 
P is the pressure of the gas and q = —nVT is the heat 
flux, with n being the heat conductivity. 

An optical birefringence owing to the planar confine¬ 
ment of diatomic molecules has been recently analyzed 
and described in terms of the expectation value of cos^ Oz , 
where Oz is the angle between the axis of a molecule 
and the propagation axis 2 : of probe light [36] . The bire¬ 
fringence signal is proportional to (l — Oz)) , be¬ 

coming non-zero when the molecules are confined to the 
xy plane Oz) < 1/3). For a gas of optically cen¬ 

trifuged molecules in the quantum state IJ, Mj = J ), 
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where Mj is the projection of J on one finds [37]: 
{cos‘^ Oz) = {J, J\ cos^ Oz\J, J) = (2J + 3)“^. As a re¬ 
sult, the birefringence signal is proportional to Since 
the rotational energy scales with the same power of J, 
one can show that the growth of the gas temperature 
due to the rotation-translation energy exchange can be 
expressed as T(r, t) = To-hAT(r) [1 — exp(—t/r^)], where 

is the decay constant of the birefringence signal, Tq is 
the ambient temperature and AT(r) is the temperature 
increase at the end of the thermalization process. The 
dependence of the latter on r reflects the gaussian profile 
of the centrifuge beam. 
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